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In a recent publication [Phys. Rev. Lett. 97, 227402 (2006), arXiv: cond-mat/0611411], it has been 
demonstrated numerically that a long-range disorder potential in semiconductor quantum wells can 
be reconstructed reliably via single-photon interferometry of spontaneously emitted light. 

In the present paper, a simplified analytical model of independent two-level systems is presented 
in order to study the reconstruction procedure in more detail. With the help of this model, the 
measured photon correlations can be calculated analytically and the influence of parameters such 
as the disorder length scale, the wavelength of the used light, or the spotsize can be investigated 
systematically. Furthermore, the relation between the proposed angle-resolved single-photon cor- 
relations and the disorder potential can be understood and the measured signal is expected to be 
closely related to the characteristic strength and length scale of the disorder. 

PACS numbers: 78.55.-m, 42.50.-p, 71.35.-y, 78.30.Ly 



I. INTRODUCTION 

Disorder in semiconductor heterostructures has 
strong influence on their optoelectronic properties. 
Independently of whether disorder is introduced by 
design or due to random interface roughness or com- 
positional fluctuations it can significantly alter the 
properties of the heterostructure [1, 2, 3, 4, 5, 6, 
]. Recent publications show that long-range spa- 
tial disorder is an important topic from both the 
technological and the fundamental point of view 
[5, 8, 9, 10]. 

In order to enhance our understanding of the role 
of disorder in these systems, we have recently pro- 
posed an experimental scheme [11] which can be 
viewed as the Fourier analogue of nano- or micro- 
luminescence [12, 13, 14, 15]. This scheme is based 
on measuring angular correlations of spontaneously 
emitted light and has been shown to give direct ac- 
cess to the spatial distribution of the optically ac- 
tive electronic states and to the effect of disorder on 
them. As our numerical investigations have shown 
the spatial distribution can be recovered very reli- 
ably via a Fourier transformation of the experimen- 
tally measurable photonic correlation [11]. 

In the present paper we present an extended ana- 
lytical study using a simplified model of uncoupled 
two- level systems. The additional simplification in- 
troduced by neglecting the many-body interactions 
allows us to deepen our understanding and present 
the fundamental principles behind the suggested re- 
construction procedure in a more transparent way. 
One of the aims is to get a more intuitive under- 
standing of the measured correlation function, such 
that a basic understanding of the characteristic pa- 



rameters of the disorder landscape can be gained 
from the direct measurement results without any re- 
construction procedure. 

The paper is structured as follows: We describe 
the experimental setup, introduce the important 
photon-correlation functions, and outline the re- 
construction procedure in Sec. II. In Sec. Ill we 
present general analytical results for a system of non- 
interacting two-level systems where the disorder ap- 
pears as a varying transition energy from site to site. 
Finally, we choose a sinusoidal model potential in- 
stead of a truly random disorder potential. For this 
special case, further analytical results are derived in 
Sec. IV before we conclude. 



II. MEASURING CORRELATIONS 

Our suggested scheme is based on the angular 
photonic correlations of spontaneously emitted light. 
Previous calculations have shown that not only pho- 
ton numbers (B^Bq), but also single-photon corre- 
lations of the form (B^Bqi) between photons of dif- 
ferent modes (q ^ q') are building up when a semi- 
conductor heterostructure reaches quasi-equilibrium 
and spontaneously emits light [16]. While the rate 
of photons ^.(B^Bq) is proportional to the photo- 
luminescence (PL) spectrum at energy ui q — c|q| 
and thus directly accessible to experiment [16, 17], 
a clever setup must be used in order to measure 
the complex correlations. The simplest possibility 
is to combine the light propagating along two dif- 
ferent directions in a common detector. Here, the 
light emitted into directions q and q' is redirected 
with the help of mirrors and through a beam splitter 
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combined into the same detector. Consequently, the 
detector detects the combined beam and the correct 
detector operator for the description of the measure- 
ment process is given by 

d q , q /=B q + e < *B q /, (1) 

i.e. a superposition of the photon operators of the 
two directions. The variable phase <f> is adjusted via 
the optical path length difference between the two 
light beams. The detected signal is proportional to 
the number operator corresponding to dq iq ' and thus 
given by 

(<, q '<W> = ((B\ + e-^B\) (B q + e^B q ,)) 
= (B q 5 q ) + (b\,B^) + 2 Re (e^ qJ B q ,) .(2) 

The first two terms on the right hand side of Eq. (2) 
are nothing but the separate PL spectra emitted into 
the two different directions while the interference 
term shows that also photon-correlations between 
the two different emission directions can be made 
visible. As previous results have shown, correlations 
can only build up if q|| = q'| in the case of perfectly 
ordered system [ ]. Here the parallel subscript de- 
notes the component of the wave vector along the 
quantum well (QW). 

While the measured signal according to Eq. (2) 
still depends on a total of six components of q and 
q', a few simplifications can be introduced accord- 
ing to the actual experiment; firstly, a spectrome- 
ter is used in order to achieve a spectrally and an- 
gularly resolved signal. The spectrometer is intro- 
duced after recombining the two emission directions 
into one beam and before the detector. Therefore, 
only the signal for equal |q| = |q'| has to be com- 
puted. Furthermore, our previous numerical inves- 
tigations have shown that the interference depends 
much stronger on the difference qy — qj| than on the 
average value such that we define the angular corre- 
lation 

U hu (M\\) = (BlB^,) (3) 

with the specific choice of q' = (qj_, — Aqy/2) and 
q = (q_L,Aqii/2) as the key observable. The mag- 
nitude of the perpendicular component of q and q' 
has to be fixed for each choice of Aqy according to 
LO 2 = c 2 (|q_L| 2 + |^Q|| | 2 /4) . Consequently, C//^(Aqy) 
depends on only two parameters, which is sufficient 
for a spectrally and angularly resolved scan. 

While the magnitude of qj^ is determined by w 
and Aqy, the sign can still be chosen differently. 
Our above choice of identical signs for the perpen- 
dicular components of q and q' represents measure- 
ment in e.g. reflection geometry as depicted in Fig. 1 
which has the advantage that non-transparent sub- 
strate can be used as both emission paths are on the 




FIG. 1: Schematic setup of a possible experimental ar- 
rangement for measuring photonic correlations. 



same side of the sample. An alternative possibility is 
to use light emitted into opposite directions which 
corresponds to q' = q. This setup was success- 
fully applied to detect angular photonic correlations 
in ordered systems [18]. 

The principle idea behind our scheme for the re- 
construction of the disorder potential exploits the 
fact that the strict momentum conservation along 
the direction parallel to the QW is broken by a disor- 
der potential. In contrast to the ordered case where 
interference could only be observed for Aqy = 0, we 
expect that non-zero correlations cannot only be ob- 
served for vanishing Aqy but also for other values. 
Please note that no such restriction applies to the 
direction perpendicular to the QW, independently 
if it is perfectly ordered or disordered, because the 
translational invariance is broken by the electronic 
confinement. 



A. Reconstruction procedure 

For the sake of simplicity of the analytical calcula- 
tions we consider a one-dimensional quantum- wire in 
the remainder of the paper. In this case, the parallel 
component qy and the spatial coordinate x along the 
wire are scalars. Apparently, all conclusions about 
momentum conservation hold as for a QW. To sim- 
plify the notation we use Ag instead of Aqy in the 
rest of the present work. During the calculations we 
limit ourselves to two emission directions with both 
planes of emission (spanned by the one-dimensional 
wire and the emission direction) being identical. 

According to Eq. (2), the correlation function 
from Eq. (3) can be experimentally measured as an 
interference contrast [11]. By varying the optical 
path length in a controlled way, it is possible to 
extract both real and imaginary part of [/^(Ag) 
such that it is legitimate to study the full complex 
value of U hlJJ (Aq). 

By Fourier transforming the measured signal ac- 
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cording to 



A 9ma3t 

~ J Ur^(A q )e- tAcix dAq 



(4) 



-A 9ma: 



one can obtain information on the spatial distribu- 
tion of emitting sources of energy Tito. Here Ag max is 
determined by the optical cone; since our scheme is 
a far-field method, the maximum value of Aq acces- 
sible in an experiment is obtained for oblique inci- 
dence and given by approximately Aq max = 2u>o jc = 
47r/Ao, where c is the speed of light in vacuum and 
Ao and tuo are the wavelength and frequency of the 
average emission frequency. 

For a fixed value of uj, Eq. (4) can exhibit multi- 
ple peaks at all those positions where the band gap 
corresponds to hu [11]. Thus, Unu,{x) can be viewed 
as a probability distribution of the local contribu- 
tions to the emission spectrum. For a fixed value of 
x, on the other hand, Un^ix) has only a single peak 
at the local transition energy such that the average 
energy can be defined as the weighted mean of the 
transition energies via 

E2 E2 

hujUhuj(x) dTiu / / Uhu>{x) ATiuj . (5) 



U(x) 



Here, E\ and Ei are the limits of the energy scan 
carried out in the experiment. U(x) is called the 
reconstructed disorder potential [11]. 

It is important to realize that the reconstruction 
procedure can only work within certain limits. Natu- 
rally, no disorder varying on length scales larger than 
the spot size can be detected. Furthermore, as with 
every far-field method, we arc limited by the wave- 
length of the emission and the reconstructed disorder 
potential can only be determined if it varies on a typ- 
ical length scale larger than the wavclenght of light. 
In the present paper we demonstrate this scheme on 
the basis of a theoretically calculated Unuj (Aq) for 
an ensemble of two-level systems. 



III. CORRELATIONS OF TWO-LEVEL 
SYSTEMS 

As starting point we consider an ensemble of un- 
coupled two-level systems as our model system. All 
two-level systems are spaced equally along x with 
lattice constant a such that the j th two-level system 
is positioned at the lattice site Rj = aj. By using 
an equation-of-motion approach and evaluating the 
Heisenberg equation of motion [11], we have shown 
previously that the correlation function can be ex- 
pressed as 



U hul (Aq) 



|2 \ ~" 



jSje iAqR i 



^( Ej -huj) 2 + ~f 2 



(6) 



Here, the prefactor contains the matrix element 
of the light-matter coupling including the micro- 
scopic intcrband dipolc matrix element, the so-called 
vacuum field amplitude, and the mode strength. 
The remaining term is a sum over sites with ho- 
mogeneously broadened peaks centered at the site 
transition energies Ej multiplied by the PL sources 
Sj = fjfj, where f? and /j l are the electron and 
hole occupation probabilities of the respective site. 
The phase factor is a consequence of the different 
emission directions and is the key ingredient for ob- 
serving the angle-resolved interference effects. 

Taking the continuum limit a — > with con- 
stant PL source density Sj/a — > S(x) and a con- 
tinuously changing transition energy E(x), we can 
rewrite Eq. (6) as 



U, 



huj 



je 



iAqx 



S(x) 



(E(x) ~ nuy + T 



dx. (7) 



In the ordered case, all transition energies E(x) — 
Eq are equal such that the angular correlation 
is proportional to a single broadened Lorentzian 
multiplied by the Fourier transform of the spot, 
/ e iAqx S(x) dx. 



A. Spatial information 

According to Eq. (4) , the Fourier transform of the 
correlation function gives some spatial information 
about the position of emitters of a certain frequency 
u. Starting with the result of Eq. (7), we obtain the 
Fourier transform 

2|^| 2 Ag max 

Ufr u (x) = 

n 



S(y) sinc(Ag max (y -x))dy (8) 



where we introduced the broadened ^-function 



1 



7 



7T E 2 + 7 2 



(9) 



and used the sine-cardinal sinc(x) = sin(x)/x. The 
broadening of g~ t is caused by the homogeneous line 
width 7 and results in a finite energy resolution. But 
also Aq max sinc(A<7 max x) approaches a Dirac delta 
function for large values of Aq max . In that case, the 
broadening of the sine function for finite values of 
Ag max reflects the fact that no spatial information 
can be obtained for variations of the site energies on 
length scales below the wavelength Ao = 47r/Ag max . 



B. Dimensionless parameters 

As one can anticipate from Eq. (8), the possibil- 
ity to obtain precise spatial information strongly de- 
pends on the ratio 7/W between homogeneous 7 and 
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the disorder strength and on the ratios \q/L and 
S/L between the observed wavelength or the spot 
size S and the typical disorder length scale L. 

In order to obtain scaled equations with respect 
to the characteristic parameters W and L of the dis- 
order potential, we introduce u = x/L and Q — 
LAq as new spatial and momentum coordinates and 
£1 = huj/W and a — 7/W as the rescaled fre- 
quency and homogeneous line width. The maxi- 
mum value of u is determined by the spotsize to be 
\u\ < u max — S /2L while Q is confined according 
to \Q\ < Q m ax = ^y ± - The dimensionless transi- 
tion energy is denoted E(x) = E(x)/W. With the 
help of these new variables, we can rewrite Eqs. (6) 
and (8) as 



Un(Q) = I g a (E(u)-tl)e^ u S(u)du, (10) 

1 Qma 



Un(u) 



7T L 



la (E(v)-fl)S(v) 
sinc(Q max (t; - u)) dv , (11) 



where we have eliminated the prcfactors by the def- 
inition U hu , = 2irL\T u \ 2 /{hW)U n and g a (E) = 
n~ l a/((E 2 +a 2 ) = Wg 7 (E). By virtue of the defini- 
tion of E, its fluctuation is on the order of unity and 
the energy resolution can be directly estimated from 
the magnitude of a. In the limit of strong disorder 
(W ^S> 7) corresponding to a — > the broadened g a 
becomes a sharp Dirac delta function and we obtain 



Un(Q) 
Usi{u) 



W»7 



W»7 



5> lQ "" S(u n ) 



dE 
du 



u — u n (O) 




1 Qn 



TV L 



/"] sinc(Q max (M n - u)) 



S(u n 



dE 
du 



(13) 



u — u n (O) 



where in both equations the sum over n runs over all 
roots u n — u n (yi) of the argument of the delta func- 
tion. These roots are defined via E(u n (Cl)) = f2. 
Obviously, only positions u n within the spot can 
contribute to the correlations, as is guaranteed by 
the factor S(u n ) which vanishes outside the excita- 
tion spot. 

Another interesting limit is obtained for a disorder 
length scale well above the wavelength of the emit- 
ted light, but still below the spotsize. In that limit 
Qmax 2> 1 such that the sine cardinal function can 
be treated as a sharp delta function and we find 



Un(u) 



1 Qr, 



£>A 



7T L 



'-g a (E(u)~n)S(u).(U) 



In that case, the spatial resolution is very good even 
for relatively large values of 7 since for fixed u = 



x/L the correlation function Ua(u) has its maximum 
at the transition energy hu = E(x). According to 
Eq. (5) the reconstruced potential follows exactly 
the actual energetic variation E{x). 



IV. MODEL POTENTIAL 

While the results of the previous section are valid 
for a periodic arrangement of non-interacting two- 
level systems with arbitrary disorder (i.e. arbitrary 
variation of the transition energies), we will intro- 
duce additional approximations in the present sec- 
tion in order to simplify the equations even further. 
In particular, we study the interpretation of Uq(Q), 
i.e., the direct measurement signal. That way, we 
can see what information on the disorder potential 
can be gained from the measurement without per- 
forming any reconstruction. 

Firstly, we assume an infinitely large spot and a 
homogeneous luminescence source S(x) = Sq. Sec- 
ondly, we consider a very specific variation of tran- 
sition energies in the form E(u) — Eq + sin(27ra). 
While a real disorder potential in general is more 
randomly fluctuating and might even contain several 
characteristic length scales, our more regular model 
potential has a perfectly well defined length scale. 

With the above mentioned approximations, we 
can simplify Eq. (10) for = Eq and get 

Un=E (Q) = [ g a (sm(2iru))e* Qu du, (15) 



After rewriting the broadened Lorentzian function 
as 



g ce (sin(27ru)) = 



1 



1 



2ty V sin(27m) + ia sin(27ru) — ia 



(16) 

we can solve the integral in Eq. (15) with the help of 
Cauchy's theorem where the contour of the integra- 
tion has to be closed in the upper (lower) complex 
plane for positive (negative) sign of Q, respectively. 

The poles of the two contributions of Eq. (16) are 
given by the roots of the energy denominators. The 
roots of the first denominator are defined by the 
equation 

sin(27ru) = — ia . (17) 
Introducing the quantity 

Ho = asinh(a) = In (a + \/l + a 2 ^j . (18) 

we can express the solutions of Eq. (17) as 



,(m) 



y a +m, u^l = +^—ya+m+^- , (19) 



2tt 



2tt' 



for integer values of m. The subscripts + and — 
signify whether the pole lies in the upper or lower 
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complex plane. Similarly, the poles of the second 
contribution of Eq. (16) are given by 



(m) I -I 



u^=+—y +m. (20) 



by 



The residues of the integrand of Eq. (15) are given 



Res 



Res 



2-7T sin(27r«) + ia ' 



2n sin(27rw) — ia ' 



= ± 



J_ 1 

2iri 27r cosh(yo) 



j_ 1 e iQui Cl 
2ni 2n cosh(yo) 



(21) 



(22) 

such that after inserting the explicit solutions for 
u i/2 ± we obtain 



Un=E {Q) 



\Q\vo 

1 e 2> 



2tt cosh(yo) 
Employing the identity 



E 



n— — og 



JQn/2 



(23) 



E 



iQn/2 



= 4tt ^ 5(Q-4ttto) (24) 



we obtain the final result 



Dh=Bb(Q) = ^ <5(Q - Anm) (25) 

9 + 9 „, „ 



Here we have introduced 

g = a + y/l 



(26) 



and used the definition of yo as well as the relations 



cosh(y ) 



(cosh(yo) 
\fl + a 2 ■■ 



■ sinh(yo)) 



9^7) 
(28) 



0.5 



cr 
< 



a) 












. J' \ 


1 


| 4n 1 L S\ 






FIG. 2: a) Correlation function Un^=E (Aq) as function 
of Aq/qo for a sinusoidal potential and different spot- 
sizes. Comparison between the analytically computed 
envelope function (thick solid line) and the numerical 
result for a spotsize of S = 10 /im (thin solid line) and 
S — 20 /im (dashed line). Furthermore indicated are 
the peakwidth proportional to the inverse spotsize and 
the peak distance inversely proportaion to the disorder 
length scale. Numerical parameters are 7 = 1 meV, 
Wh = 20, L = 1 mm 

b) The height of the second maximum of Un^=E {A.q) 
as a function of the variance of the disorder potential. 
Thick line: sinusoidal potential, analytical result. Dots: 
numerical results for arbitrary spot size. The dashed line 
gives the thick line divided by 2. For random potentials 
of various length scales L the resulting curves are situ- 
ated in the shaded area. The thin line is a guide for the 
eye. 



Our result, Eq. (25), nicely demonstrates that for 
the sinusoidal model potential the expected angular 
correlation consists of regularly spaced delta func- 
tions multiplied with an envelop function which ex- 
ponentially decays with \Q\. The decay rate with 
respect to \Q\ is determined by ln(g) and thus a func- 
tion of the ratio a — j/W. The different peaks are 
located at integer multiples of 4ir corresponding to 
a spacing of 47r/L for the unsealed variable Aq. 

In order to confirm our results and study the in- 
fluence of a finite spot size, we display the envelope 
function 4g~^ /(g + g^ 1 ) together with the results 
of two numerical computations for finite spot size 
in Fig. 2a. Here, we have normalized the result to 
a maximum value of 1 and plotted the correlation 
as a function of Aq/q where the wavenumber qg 
corresponds to an assumed emission wavelength of 



A = 800 nm. The distance between neighbouring 
sites has been chosen a = 5 nm. 

We notice that our predictions are confirmed by 
the numerical simulations and that furthermore the 
finite spotsize does not influence the envelope func- 
tion. While the width of the peaks is inversely pro- 
portional to the spotsize, both their height and their 
position do not depend on S. In a real experiment, 
we can thus expect that the distance between the 
central peak and the first neighbouring peak is a sig- 
nature of the longest disorder length scale involved, 
while the ratio between the peak heights of first and 
second peak, given by 



<7n=,E; (47r) 
Ua= Eo (0) 



= 9 



(29) 



G 



can provide a measure of the ratio between homoge- 
neous 7 and the disorder fluctuation W. For known 
homogeneous broadening, one can thus extract in- 
formation on the energetic spread of the transition 
energies. 

A good measure of the energy fluctuations for very 
different kinds of disorder potentials is given by the 
standard deviation a. For a general random poten- 
tial used in a numerical computation, the standard 
deviation is obtained via 

3 

as a sum over the lattice sites, where E is the average 
energy and N is the total number of sites. For our 
sinusoidal model potential in the continuum limit, 
we must use the integral form 

o 2 = \ \ sin 2 (27ra;/i) dx = W 2 /2 . (31) 
L Jo 

The result shows that the standard deviation in 
that case is simply proportional to the amplitude 
of the sinusoidal function. For a comparison be- 
tween the analytical model system and numerical 
solutions with random disorder potentials, a sys- 
tematic study of the ratio g~ 2 as function of cr/7 is 
shown in Fig. 2b. As a consequence of the random 
disorder potential, also the peak ratios for different 
realizations are different. The shaded area denotes 
the region in which all numerical results of a large 
number of simulations have been found. The ana- 
lytical result (thick solid line) is found to be roughly 
twice as large over the whole range of a values. An- 
alytical tests with a superposition of more than one 
sine function have shown that the analytical result 
for such a model also decreases towards the numer- 
ically observed region. 



V. CONCLUSIONS 

In conclusion, the presented analytical derivation 
of angular photonic correlations supports the results 
of previously published numerical investigations [11] 
and helps to understand the underlying principles of 
the suggested reconstruction scheme. 

We have derived the angular correlation function 
for an ensemble of uncoupled two-level systems. 
In that model, the correlation function has been 
found to be given by the Fourier transform of the 
excitation spot for a perfectly ordered system while 
the presence of disorder shows up in a deviation 
from this shape. The spatial distribution of emitters 
can be recovered via a Fourier transformation. For a 
sinusoidal model potential analytical formulas have 
been derived for the dependency of the primary ex- 
perimentally detectable signal on the characteristic 
strength and length scale of the disorder potential. 
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